Robust plasmon waveguides in strongly-interacting nanowire arrays 
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Arrays of parallel metallic nanowires are shown to provide a tunable, robust, and versatile platform 
for plasmon interconnects, including high- curvature turns with minimum signal loss. The proposed 
guiding mechanism relies on gap plasmons existing in the region between adjacent nanowires of 
dimers and mult i- wire arrays. We focus on square and circular silver nanowires in silica, for which 
excellent agreement between both boundary element method and multiple multipolar expansion 
calculations is obtained. Our work provides the tools for designing plasmon-based interconnects 
and achieving high degree of integration with minimum cross talk between adjacent plasmon guides. 
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Electromagnetic modes in metal surfaces known as 
plasmons can propagate along millimeters in metallic 
structures at near-infrared frequencies [1], thus providing 
a plausible substitute for the electrical impulses used in 
current electronic circuits operating at microwave clock 
frequencies 0, Q. Several designs of plasmon intercon- 
nects have been prototyped in recent years, including 
metallic waveguides of finite cross section in symmet- 
ric 0, [H[ and asymmetric [6] environments, channels cut 
into flat surfaces [7[ , plasmon-band-gap structures based 
upon periodic corrugations |8], and plasmon hopping in 
arrays of nanoparticles 0, Plasmon modes can be 
tuned in frequency, and their spatial distribution molded, 
by tailoring the geometry of metallic structures on the 
nanometer scale. In particular, extreme plasmon confine- 
ment has been achieved in narrow insulator films buried 
inside metal [ll|. Actually, buried structures provide a 
natural but technologically challenging approach to com- 
pact integration. In contrast to that, open plasmonic 
geometries involve electromagnetic fields extending sig- 
nificantly away from the metal @, @, 0, H, 0, E3], and 
consequently producing a substantial degree of cross-talk 
between neighboring waveguides [12[ . 

In this Letter, arrays of parallel metallic nanowires 
are shown to provide a versatile and tunable platform 
for highly-integrated plasmon interconnects. The prop- 
agation distance and degree of confinement of the plas- 
mon guided modes depend strongly on the separation 
between wires. Individual wire modes are recovered at 
large separations, while mode hybridization is observed 
when the spacing is reduced. Gap modes are observed at 
small separations, highly localized in the regions between 
two adjacent wires. We use both the boundary element 
method (BEM) [13] and a two-dimensional multiple- 
elastic-scattering multipolar expansion of the fields for 
straight cylinders (2D-MESME) Q, with the two ap- 
proaches resulting in complete agreement on the scale 
of the plots. These methods provide rigorous solutions 
of Maxwell's equations in frequency space for materials 
described by local dielectric functions and separated by 
abrupt interfaces. We focus on silver nanowires of cir- 



cular and square cross sections embedded in silica. The 
dielectric functions of silver [15[ and silica [16[ have been 
taken from optical data. The proposed guiding mech- 
anism is demonstrated to be tolerant to asymmetry in 
wire dimers and sharp turns of subwavelength radius. 

The localized plasmons sustained by our structures can 
be conveniently characterized using the photonic local 
density of states (LDOS), defined by analogy to its elec- 
tronic counterpart as the combined local intensity of all 
eigenmodes of the system under investigation. We in par- 
ticular consider the LDOS relative to its value in vacuum, 
uj 2 /3tt 2 c 3 . The LDOS is proportional to the radiative de- 
cay rate of excited atoms [17j], which we in turn obtain 
using BEM from the imaginary part of the self-induced 
electric field acting back on a dipole [18[. We have dou- 
ble checked our results by comparing with the excess of 
space-integrated total density of states (DOS) with re- 
spect to vacuum, which is directly accessible through 2D- 
MESME 0. 

We start by considering a dimer formed by two 200- 
nm circular silver wires embedded in silica, as shown in 
Fig. QJa). The contour plots of Fig. QJb)-(e) show the 
DOS resolved in contributions of different momentum fcy 
parallel to the wires for various separations between wire 
surfaces (d). A strongly bound mode is observed at small 
wire separations [Fig. [2(b)] , with fey well above fch? the 
momentum of light in the host silica. The spatial ex- 
tension of this gap mode is limited to the inter-wire re- 
gion (see inset), and thus, it is expected to interact very 
weakly with other structures sitting in the vicinity of 
the wires but far from the gap. This mode evolves con- 
tinuously for increasing inter-wire distance to become a 
hybridized monopole-monopole mode of induced-charge 
pattern (+) •••(—) aligned with the dimer axis [Fig. 
njc)-(d)]. This is in contrast to the (H — ) • • • (H — ) 
dipole-dipole plasmon in particle dimers [19], which is 
the lowest-energy mode according to plasmon chemistry 
arguments [20[. In this sense, wires are distinctly differ- 
ent from particles because charge neutrality is guaranteed 
by oscillations along the rods for finite fcy, thus making 
two-dimensional monopoles possible. At sufficiently large 
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FIG. 1: Gap plasmon modes of two parallel silver nanowires 
in silica, (a) Schematic view of the geometry, (b)-(e) Pho- 
tonic density of states (DOS) as a function of energy and mo- 
mentum parallel to the wires for various dimer separations d. 
The insets show the spatial distribution of the local density of 
states (LDOS) for the lowest-energy gap mode at a free-space 
light wavelength of 1550 nm. Brighter regions correspond to 
higher DOS and LDOS. The maximum LDOS in the inset of 
(c) is ~ 15000 times the vacuum value. 



SiCX 




A Ag M Ag ■ Ag 1 Ag 



10 nm 



200 nm 



f (b) 



> 
s— 

CD 
C 

CD , 

c 
o 

o 




10 



20 



30 




Parallel momentum (jam ) 



FIG. 2: Evolution of gap plasmon modes with the number of 
nanowires in an array (N). (a) Schematic view of the geome- 
try. (b)-(e) DOS as a function of energy and momentum k\\ 
parallel to the wires for arrays of N = 1,3, and oo nanowires. 
The insets show LDOS maps for the lowest-frequency gap 
mode at a wavelength of 1550 nm. The lower inset in (e) 
shows the transverse momentum dependence of the DOS for 
k\\ = 16 /im" 1 compared to a tight-binding model of the gap 
mode (solid curve) over the first Brillouin zone of the ID lat- 
tice. 



distance, single-wire plasmons of m = azimuthal sym- 
metry are recovered [cf. Fig. [TJe) and Fig. |2jb); see Ref. 
[2lj for analytical expressions of single- wire plasmons]. 

Fig. [2] proves that the gap mode is really local. The 
lowest-frequency plasmon branch of the trimers in Figs. 
[21(c) and [2^d) follows approximately the dispersion rela- 
tion of the gap mode in the dimer with the same gap 
distance d — 10 nm [Fig. HJc)], although closer exami- 
nation reveals the splitting of this mode into two very 
close modes [see lower inset in Fig. Efd)]. In the infinite 
wire array of Fig. [2][e) a plasmon branch of gap modes 
is found for each value of the transverse momentum k t , 
parallel to the array plane and perpendicular to the wires 
(k t = in the figure). Incidentally, propagation across 
wires mimics plasmon hopping in particle chains [221 ]. 

Nearly-touching trimers can be regarded as dimers 
formed by two coupled gaps, and similarly, AT- wire arrays 
behave as structures formed by N — 1 gaps. The large 
degree of plasmon localization observed for small d sug- 
gests using a tight-binding model [23[ , in which an unper- 
turbed Hamiltonian Ho describes uncoupled gap modes 
\j) at sites j, such that (J\H \f) = u k where cu k 



is the mode energy for fixed fen. Neighboring gaps can 
interact in this model via a potential V, with non-zero 
matrix elements (j\V\j ± 1) = where A^ is the 

hopping energy. Then, the two linear-trimer modes have 
energies lj^ ± A^/2 [i.e., the dimer band lies halfway 
between the two trimer bands, which we have verified 
by comparison of Figs. [TJc) and EJd)]. Also, the plas- 
mon modes of an infinite, periodic wire array must have 
the form \^k t ) = Xlj ex P(i^£ a i) \ j) m virtue of Bloch's 
theorem [23], where a is the period. These states diag- 
onalize the full Hamiltonian Ho + V and have energies 
LJk\\k t ~ + Afen cos(k t a). We have tested this formula 
in the array of Fig. 2(e), where the lower inset shows 
LJk\\k t as a function of k t for fey = 16/im _1 (solid curve) 
compared with the actual 2D-MESME calculation of the 
DOS. The unperturbed gap energy lj^ = 1.18 eV and the 
hopping parameter A^ = 0.01 eV have been taken from 
the dimer [Figs. HJc)] and trimer [Figs. EJd)] with the 
same value of d, respectively. Similar agreement between 
model and full calculation is observed over the range of 
fen under consideration. 
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FIG. 3: (a) Phase velocity of gap plasmon modes in silver 
wire dimers and infinite- arrays (for kt = and kt — tv/cl) as 
a function of separation d for fixed wavelength A = 1550 nm. 
(b) Propagation length L under the same conditions as in (a), 
obtained from L = l/2Im{/c||}, where Im{/c||} corresponds to 
the HWHM of the ku -dependent DOS. 



Interestingly, these analytical expressions apply to the 
large d limit as well, in which the tight-binding model is 
constructed based upon localized plasmons of the wires, 
showing similar agreement with rigorous DOS calcula- 
tions. Tight-binding is thus the natural description of 
both the small and large d limits in the noted continu- 
ous evolution from the localized gap mode (small d) to 
the lowest-energy hybridized (monopole-monopole) wire 
modes (large d) [19[- 

The degree of plasmon localization increases with de- 
creasing gap distance d. This is reflected in a reduction 
of the phase velocity v p = Chkh/k\\ relative to the speed 
of light in silica, c^, as shown in Fig. [31(a). The group 
velocity (no shown) is also reduced, so that gap modes be- 
come considerably slower than light in silica. The prop- 
agation distance is strongly-dependent on inter-wire dis- 
tance [Fig. [31(b)]: the large confinement observed at small 
separations increases the relative weight of the electric 
field intensity inside the metal, where ohmic losses are 
produced in proportion to that intensity within linear 
response. Nevertheless, the gap mode involves electric 
field polarization mainly perpendicular to the wire sur- 
faces near the gap, where light energy is concentrated, 
and this is beneficial to obtain longer propagation dis- 
tances because the normal electric field inside the metal 
is reduced by its large dielectric function to fulfill the con- 
tinuity of the normal electric displacement. This gives 
rise to propagation distances of the order of tens of mi- 
crons for separations of tens of nanometers, accompanied 



by relatively large mode confinement. The tradeoff be- 
tween confinement and propagation distance is clearly 
illustrated in the long-<i behavior of the infinite array for 
kt = and k t = n/a, with the former showing longer 
propagation and smaller phase velocity (see Fig. [3]). 
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FIG. 4: Gap mode in co-planar and co-axial bi-tori compared 
with a straight- wire dimer for a gap distance d — 10 nm. Par- 
tial contributions to the LDOS are shown as a function of 
wavelength for a point in the center of the gap, both at fixed 
azimuthal number in tori (ra = 8 for radius b — 750 nm and 
m — 4 for b = 375 nm) and at fixed parallel momentum in the 
straight dimer (fey ~ 10.7 /mi -1 , such that k\\ — m/b). The 
curves are normalized to their maximum value. 

We find it convenient to define a figure of merit F for 
the waveguides expressed as the ratio between the prop- 
agation distance and the geometric mean of the mode 
diameter in the transverse directions. The quantities 
F 2 and F 3 should be roughly proportional to the num- 
ber of logical elements that can be integrated using a 
given waveguiding scheme with two and three dimen- 
sional packing, respectively. We obtain F ~ 540 for the 
wire dimer of Fig. [3] at a separation of 10 nm. This has to 
be compared with values of F < 50 for channel plasmon 
polaritons [7[ and particle arrays [10[. We conclude that 
wire arrays yield high values of F, also improved with re- 
spect to those obtained for single wires (e.g., F « 100 at 
100 nm radius and 1550 nm wavelength). The decrease in 
propagation length is the price to pay for plasmon con- 
finement, but wire arrays seem to perform optimally with 
respect to the figure of merit F. 

Reliable plasmon waveguides must be robust against 
fabrication imperfections and sharp turns. Next, we show 
that gap waveguides satisfy these requirements. In par- 
ticular, curved waveguide paths produce radiative losses 
originating in coupling to propagating light waves when 
translational invariance is broken. We analyze this effect 
in Fig. [H both for non-identical co-planar tori and for 
identical co-axial tori, using the prescription fcy = m/b 
to compare with straight waveguide modes, where b is the 
toroidal radius (see insets) and m is the azimuthal mo- 
mentum number. The calculations are performed using 
BEM, specialized for axially-symmetric geometries [13|- 
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Radiative losses are still small compared to absorption for 
b = 750 nm (cf. curves for straight wires and large-radius 
tori in Fig. 0J showing only ~ 5% increase in peak width 
of curved versus straight wires due to radiative losses in 
the former), but they become sizable for shaper turns 
(the width increases by 40% and 95% for b = 375 nm in 
co-axial and co-planar torii, respectively). 
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FIG. 5: Gap mode against variations of wire radius (a) and 
shape (from circular to square cross section) (b). The LDOS 
is represented as a function of k\\ for a point at the center of 
each dimer and a wavelength A = 1550 nm. One of the wires 
in the dimers of (a) has a fixed radius of 100 nm, while various 
values of the radius are considered for the neighboring wire: 
100 nm, 90 nm, and 80 nm, from top to bottom. The distance 
between wires is d = 10 nm in all cases. The horizontal diam- 
eter of the wires in (b) is 200 nm for all cross sections. 

Guided gap plasmons are also robust against dimer 
asymmetries, as shown in Fig. [5fa) for fixed wavelength 
A = 1550 nm and gap distance d = 10 nm. Variations 
of up to 20% in the relative radius of neighboring wires 
produce just a small, tolerable shift in fey. However, wire 
shape is a critical parameter, which we study in Fig.[5jb) 
through the transition from circular to square cross sec- 
tion. This produces a shift of the gap plasmon towards 
larger fey, consistent with the higher degree of confine- 
ment that occurs when evolving from the line-like con- 
tact of the circular wires to the planar waveguide de- 
fined by the square wires, the guided modes of which 



have been the subject of recent experimental investiga- 
tion [11]. This increase in confinement is accompanied 
by peak broadening originating in larger overlap of the 
gap mode with the metal (ohmic losses). The observed 
extreme sensitivity to shape and separation of the wires 
imposes severe limits to the precision required in the fab- 
rication of the arrays in order to maintain a homogenous 
mode wavelength along the waveguide. 

In conclusion, we have shown that gap plasmon 
modes existing in the region defined by two neighbor- 
ing nanowires are excellent candidates to guide signals 
over tens of microns. These modes are quite robust 
against both unintended variations of wire cross section 
and curvature in short turns, and thus, gap plasmons can 
be guided with minimum losses over complicated wind- 
ing paths of micrometer dimensions. Furthermore, gap 
modes are highly confined to the gap region, so that inter- 
mixing between neighboring wire-dimers can be mini- 
mized, thus preventing waveguide cross-talk and allowing 
highly-integrated plasmonic circuits in three dimensional 
spaces. 
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